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'jrt ! In 1877, E. I. Zolotarev ^SJ |2j found an explicit expression, in terms of 

elliptic functions, of the rational function of given degree m which is uni- 
formly closest to sgn [x) on the union of two intervals [—1, —a] U [a, 1]. This 
result was subject to many generalizations, and it has applications in electric 
^ ! engineering. 

^ \ Surprisingly, to the best of our knowledge, the similar problem for poly- 

cn I nomials was not solved yet, so we investigate it in this paper. 

Q ■ For comparison, we mention here the results on the uniform approxima- 

\0 • tion of |x|", a > on [—1,1]. Polynomial approximation was studied by 

^ ! S. Bernstein |3J |3] who found that for the error Em{oi) of the best approxi- 

^-^ I mation by polynomials of degree m the following limit exists: 

c^ 

lim m°'Em{a) = ^{a) > 0. 



X 



c^ 



This result for a = 1 was obtained by Bernstein in 1914, and he asked the 



;_i ■ question, whether one can express /i(l) in terms of some known transcenden- 

tal functions. This question is still open. Bernstein also obtained in |4| the 
asymptotic relation 

limu(a) = l/2. 

a— >0 

The analogous problem of uniform rational approximation of |a;|", a > 
on [—1, 1] was recently solved by H. Stahl jl^, who completed a long line of 
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development with a remarkably explicit answer: 

lim exp(n^/arn)EL, = 2^^° I sin(7ra/2)|, 

where i?^ is the error of the best rational approximation. 

Now we state our results. Let Pm be the polynomial of degree at most 
2m + 1 of least deviation from sgn (x) on X{a) = [— 1, — a] U [a, 1] where 
< a < 1. It follows from the general theory of Chebyshev that such 
polynomial is unique. Put Lm{a) = maxx(a) \Pm.{x) — sgn (x)|. Then we have 

Theorem 1 The following limit exists 

l + aX"" , , , 1-a 



lim ^/m -„.y-~, I — 



Remark. When approximating an odd function on a symmetric set, poly- 
nomials of even degrees are useless. Indeed, if g is a polynomial of degree at 
most 2m which deviates least from our function, then (^(2;) — q{—z))/2 is an 
odd polynomial, thus its degree is less than 2m and its deviation is at most 
that of q. Thus q is of odd degree. 

Our approximation problem is equivalent to a problem of weighted ap- 
proximation on a single interval. Indeed, Pm can be written as Pm{x) = 
xqm{x^)i where g^ is the polynomial of degree m that minimizes the weighted 
uniform distance 

sup y/x\q{x) - l/\/x\. (1) 

over all polynomials q of degree at most m. 

It is useful to compare our result with the result of Bernstein [0] , see also 
[H Additions and Problems, 44] that gives the rate of the best unweighted 
uniform polynomial approximation of 1/ \/x: 

lim Jm { I inf sup \r{x) — 1/ \fx\ = — ^(1 — a^)a~^'^. (2) 

m^oo^ \l-aj degr=m [„2^^] ' ^ ' '^ ^ 2^^ ' ^ ' 

In our proof of Theorem 1, the asymptotics of the error term is obtained 
in the form 

„„,Vis(i±^)'"L„W = e-^Sl^, (3) 

'^^'^ \\- a) xfa 



with 

where X[2,oo) is the characteristic function of the ray [2, oo), and H is the 
conformal map of the upper half-plane onto the region in the upper half 
plane above the curve 

{t + i arccos e^* : t > 0}, 

normalized by -ff(O) = and H{z) ^ z^ z ^ oo. Then the numerical value 
c = (1/2) log(27r) is derived from comparison of Q with Q for a — > 1. So, as 
a curious corollary from Theorem 1 and the result of Bernstein, we evaluate 
the integral (JD). 

Following Bernstein, we also consider approximation by entire functions 
of exponential type. Let i^(A) be the error of the best uniform approximation 
of sgn (x) by entire functions of exponential type one, on the set (— cxd, —A] U 
[A,+oo). 

Theorem 2 The following limit exists 

lim ^/Aexp{A)L{A) = \Fl/^. 



The proof of Theorem 2 is similar to (and simpler than) that of Theorem 1. 
On the best L} approximation of sgn (x) by entire functions of exponential 
type we refer to [16]. 

Our proofs are based on special representations of polynomials and entire 
functions of best approximation which are of independent interest. It follows 
from Chebyshev's theory (see, for example, pi Ch. II]) that polynomials pm 
are characterized by the property that the difference Pm{x) — sgn (x) takes 
its extreme values ±Lm(a) on X{a) 2m + 4 times intermittently, so that the 
graph of Pm looks like this: 



-1 -0.5 





Fig. 1. Graph of p4 with a = 0.1. 

The extremal entire function is unique and is characterized by the prop- 
erties that it has no asymptotic values, all its critical values are real; those 
on the negative ray are —1±L{A) and those on the positive ray are 1±L{A). 

We mention a general theorem of Maclane ^l] and Vinberg [iBi on the 
existence and uniqueness of real polynomials and entire functions with pre- 
scribed (ordered!) sequences of critical values. In the case of polynomials, 
this theorem says that there is a one to one correspondence between finite 
"up-down" real sequences 



< Cfc-l > Cfc < Cfe+i > 



and real polynomials whose all critical points are real, modulo a change of 
the independent variable z \-^ az + b, a > 0, 6 G R. If (xk) is the sequence of 
critical points of such polynomial, then Ck = P{xk). 

The MacLane-Vinberg theorem is based on an explicit description of 
the Riemann surfaces spread over the plane of the inverse functions P~^. 
This explicit description and our previous work [7j, [T3] suggest to look for 
a representation of these extremal polynomials and entire functions in the 
form cos (f){z) where is an appropriate conformal map. 

Let L and B be positive numbers, < L < 1, and 



chB 



2e" 



B 



oo. 



(5) 



Consider the component 7^ C {-2 : ^z G [0,7r], '^z > 0} of the preimage of 
the ray 

{w:^w = 1/L, 3w < 0} 

under w = cosz. It is easy to see that this curve 7^ can be parametrized as 

1B = {arccos(chi?/cht) + it : B <t < 00}. 

This curve begins at iB and then goes to infinity approaching the hne {it/2 + 
it, t > 0} with exponential rate. 

Let Qb be the region in the upper half-plane whose boundary consists 
of the positive ray, the vertical segment [0,i-B] and the curve 75. For fixed 
-B > 0, let (pB be the conformal map of the first quadrant onto Qb such that 
(j)Biz) ^ z, z -^ 00 and 0^(0) = iB. Let A = A{B) = 0;^^(O). Then A is a 
continuous strictly increasing function of B, and we may consider the inverse 
function B{A). 

Theorem 3 The approximation error in Theorem 2 is L{A) = 1/ch B{A), 
and the extremal function can be defined in the first quadrant by the formula 

1 - L{A) cos (f)B{A)- 



Let fiB,m. be the region in the half-strip {z : ^z E (0, 7r(m + 1)), '^z > 0} 
bounded on the left by 75. Let (f)B,m be the conformal map of the first 
quadrant onto i^B,m such that 0B,m(O) = iB, 0B,m(l) = 7r(m + 1) and 
0B,m(oo) = 00. Let a = a{B,m) = 0]j\„(O). Then a[B,m) is a continuous 
increasing function of B for fixed m, so it has the inverse Bm{a). 

Theorem 4 The error term in Theorem 1 is Lm{a) = l/ch Bm{a), and the 
extremal polynomial is given in the first quadrant by 

1 — Lm(a) cos (pBm, whcrc B = Bm(a). 



Discontinuous functions cannot be uniformly approximated by polynomi- 
als with arbitrarily small precision, however, in our situation we can obtain 
an approximation which seems to be the second best thing to the uniform 
approximation. 



We introduce the notation^ 

L(/, g) = mi{h : f{x - h) - h < g{x) < f{x + h) + h, -1 < x < 1}. 

Thus the statement that L(sgn ,g)<e means that the graph of the restriction 
of gf on [—1, 1] belongs to a "rectangular corridor" of width e around the 
"completed graph" of sgn {x), which consists of the graph of sgn (x) and the 
vertical segment [(0, —1), (0, 1)]. 



Fig. 2. Levy's neighborhood of the function sgn . 

It is easy to see that if a = Lm{ci) then our polynomial pm from Theorem 1 
is the unique polynomial of degree 2m + 1 which minimizes L(sgn,p). We 
have 



Theorem 5 



,. m 
lim 

m^oo log m 



L(sgn,p^ 



1 
2' 



Remarks. One could also use the Hausdorff distance between the com- 
pleted graphs. In the case of sgn (x), the Hausdorff distance will differ from 
L(sgn,p,„) by a factor of \/2. Approximation of functions with respect to 
Hausdorff distance between their completed graphs was much studied by 
Sendov [12] and his followers. An arbitrary bounded function on [0, 1] can be 
approximated in this sense by polynomials of degree n with error 0((log n)/n) 



^If [—1, 1] is replaced by (— oo, cx)) this becomes the Levy distance. It is really a distance 
on the set of bounded increasing functions on the real line 10, Ch. VIII]. 



Proof of Theorems 3 and 4- Let be either 0^ or (f)B,m- By inspection 
of the boundary correspondence, we conclude that / = 1 — Lcos0 is real 
on the positive ray and pure imaginary on the positive imaginary ray. So 
/ extends to an entire function by two reflections. The extended function 
evidently satisfies 



f{z) = f{z) and -f{-z) = f{z), 

so we conclude that / is odd. In the case of Theorem 1, the region QB,m 
is close to the strip {z : ^z E (7r/2,7r(m + 1))} as '^z — > cxd, so 4>B,m ~ 
i{2m + l) log z, z ^ oo, so / is a polynomial of degree 2m + 1. In Theorem 2, 
(f){z) ^ z, so f has exponential type one. In both theorems, differentiation 
shows that the only critical points of / in the closed right half-plane are 
preimages of the critical points of the cosine under 0. So the graph of / 
has the required shape. That a polynomial with such graph is the unique 
extremal for Theorem 1 follows from the general theorem of Chebyshev on 
the uniform approximation of continuous functions ^1 Ch. II]. 

The proof that the entire function / we just constructed is the unique 
extremal for Theorem 2 might not be so well-known, so we include this proof 
which we learned from B. Ya. Levin (compare [71ll5j). 

All critical points of our entire function / are real. Let xi < X2 < . . . be 
the sequence of positive critical points of /, Then we have xi > a, and 

/(Xfc) = 1 + (-l)'=-iL, and fiA) = l-L. (6) 

Let g be another real entire function of exponential type 1 such that 

sup \g{x) — 1| < L for X > A. (7) 

We may assume without loss of generality that g is odd (otherwise replace it 
by {g{x)—g{—x))/2 which also satisfies ((Tj)). Equations (0) and ((7j) imply that 
the graph of g intersects the graph of / on every interval [xk,Xk+i], k > 1. 
More precisely, there is a sequence of zeros y^ oi f — g (where multiple zeros 
are repeated according to their multiplicity), which is interlacent with Xk, 
that is 

Xl < yi < X2 < 1/2 < ■ ■ ■ , 

and in addition to those pj, f — g has at least one zero in (0, xi]. By the well- 
known theorem jOj VII, Thm. 1], it follows that the meromorphic function 

1 -j-j- 1 - z/xk 

kez\{o} 



^w-.,.n^,r_./,. 



has imaginary part of constant sign in the upper half-plane, and of opposite 
sign in the lower half-plane. This implies that 

F(re'^) = 0{r), (8) 

when r -^ +00, uniformly with respect to 6 for e < 6 < tt — e, for every 
e > 0. Similar estimate holds in the lower half-plane. 

As (/ - g){yk) = 0, we have 

^ = P/F, (9) 

where P is an entire function of exponential type. 

It is easy to see that the left hand side of © is bounded for \'^z\ > 1. 
Indeed, Phragmen and Lindelof give \f{z) — g{z)\ < Ci exp \'^z\, while /' has 
only real zeros and approaches Lcos(x ± a)) as x ^ ±00, where a is some 
real constant. It follows that |/'(2;)| > C2exp |Sz| for \Qz\ > 1, so (/ — g)/f' 
is bounded for \'^z\ > 1. 

So we conclude from (jHJ that P{z) = 0{\z\) and this contradicts the fact 
that P has at least two zeros, unless P = f — g = 0. This completes the 
proof. 

Proof of Theorem 1. We recall that Bm = Bm{a), ^m = ^s„,m and 
the conformal maps of the first quadrant Q onto Qm were defined before 
Theorem 4. We are going to prove Q first, which is the same as 

Bm= {m + -1 log ^— ^ + - logm + - log ^ _ ^^ + c + o(l), (10) 

as m ^ 00 and a is fixed. Here c is an absolute constant. 

We need some auxiliary conformal maps. Let tp : Q ^ Q he defined by 

^(z) = 

where 

Am= {m+-) log- >0 (11) 

The reason for such choice of Am will be seen later. Then ip gives the following 
boundary points correspondence: 

ip : (0, a, 1, cx)) h-> {iAm, 0, 00, iAm/a). 




The function $m(^) = "^^m ° ip~'^{—iz) maps the second quadrant onto i^mi 
and sends the positive imaginary axis onto the interval I = {0,iTT{m + 1)). 
By reflection, we extend $m to a map from the upper half-plane onto the 
region i{flrn U Qm U i). This extended map $m gives the following boundary 
correspondence: 

where we set Cm = ^m/o- 

Now we introduce the map hm{z) = ^m{z + ^m) — An- Our first goal 
is to show that the sequence h^ tends to a limit, and to describe this limit. 
The boundary correspondence under hm is this: 

We represent hm as the Schwarz integral of its imaginary part: 

hm{z) = U+ - log- ^-^ (-^ + - / -)vmit)dt, 

\ 2/ 1 — za/{{l — ajAmj nJ-ooKt — z tj 



where 



(12) 

I ^/^; t T- [~^m " Ami ^m ~ ^mj • 



Our choice of Am in ()lip implies that the first summand in the right hand 
side of (fT^ has a limit 

/ 1\ l + za/{{l + a)Am) 

hm m H — log — - = az, 

— V 2) '^l-za/{{l-a)Am) 

where 

^ = (l-a2)log((l + a)/(l-a))- ^^^^ 

It is easy to see that the integral in (fT^ converges to a bounded function of 
the form 

where p^- is a bounded positive function which we will describe shortly. The 
image of hm has a limit Vt* in the sense of Caratheodory; Vt* is the region in 
the upper half-plane above the graph of the function arccose"^', x > 0. So 
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hm -^ H„, where if^ is the conformal map of the upper half-plane onto ^2*, 
-f^o-(O) = and H„{z) ~ az as z — > oo. 
Thus we have a Schwarz representation 

H^{z) = az + - f (- -) pait)dt, 

Ti J \t — Z t J 

and Pcrif) = '^Hcrit) for real t. Our next goal is to study asymptotics of 
Brn = —hm{—Am). We use the following comparison function 

g^{z) = U+ - log- '-^ f-^ + - / -)xrnit)dt, 

\ 2/ 1 — za/[{l — a)Am) 2 Jo \t — z tj 

where Xm is the characteristic function of the set (— oo, — 2(Am + l)]U[2, +oo). 
We have 

1 f°° ( vr \ dt 

lim (/l„(-A^) - gm{-Ara)) = / Pait) - -X[2.oo)(t) — , 

m^oo vr JO V Z J t 

and using (fTT|l . 

gm{-Am) = -Arn - - \og{Am + 1). 

Combining these two equations we obtain 

1 1 f°° f vr \ dt 

Bm = -h^m{-Am) = Am + -log A^ + - J {Pa{t) " ■^X[2,oc>){t) j — + o(l). 

using the evident transformation law Hc,{\z) = H„\{z) we obtain Pa{Xt) = 
Pax{t), and therefore, 

1 11 f°° f TT \ dt 

5^ = A^ + - log A™ + - logo- + - / pi(t) - -X\2.oo){f) — + o(l). 

/ / TTJOV I It 

Substituting the values of Am and a from (fTTj) and (fT^ we obtain (fTUI) with 

c=-l (piW - -X[.oo)j y. (14) 

The numerical value 

c=(l/2)log(27r) 

is obtained from comparison with Bernstein's result ^. 
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Indeed, in view of (|T]), (0) and (fTTH) we have 

Lm{ci) = inf sup \^\q{x) — l/^/x\ ~ 2e^^" 

deg q=m [„2^i] 

, ^ , ,, 1 — a /I — a 
= 2(e-^ + o(l)) 



^/2a \1 + aj ^/m' 

Comparing this expression for Lm{a) with the expression Q when a ^ 1, 
we obtain ^ with c = (1/2) log(27r). 

Proof of Theorem 2. We have to prove that B = A+{l/2) log A + co + o(l) 
as A — >• oo, where Cq = (1/2) log(7r/2). 

Let /i (2;) = y/z'^ — A^ be the conformal map of the first quadrant Q onto 
itself, sending A to 0, and f\{z) ~ 2 as 2 -^ 00. Then /i(0) = iA. Let /2 be 
the conformal map of Q onto fi^, /2(0) = and /2(-2) ~ 2;, as 2; ^ 00. We 
extend /2 by symmetry, reflecting both domains in the positive ray. So from 
now on /2 is defined in the right half-plane. The condition that f2{iA) = iB 
defines the number B uniquely. It is easy to see that 

B ^ A 

as y4 ^ 00. 

Now put h{z) = if2{—iz) for convenience. This h maps the upper half- 
plane onto a subregion of the upper halfplane. The boundary of this subre- 
gion is asymptotic to the line {'^z = 7r/2} as \^z\ — > 00. 

We have B = h{A), and we wish to find asymptotics of h{A) as A — »• 00. 

To do this, we use the following comparison function 



00 



The integral in the right hand side is the Schwarz formula for an analytic 
function in the upper half-plane whose imaginary part equals 

(^/2)X(-oo-A-2]U[A+2,oo)- 

Our function h has a similar representation in terms of its imaginary part 
on the real line. Subtracting these two representations, we obtain 

2A /•«= v{t + A) 

(JjL ^ 



g{A) - h{A) = — / 
n Jo 



2At + 12 
11 



where v{t) = Q{g{t) — h{t)). Now we claim that v{t + A) tends to a hmit 
Vq, as a ^ +00. This hmit is (7r/2)x[2,oo)(^) — '^H(t), where and H{z) = 
liinA^+oo{h{z + A) — B{A)) is the conformal map of the upper half-plane onto 
the region in the upper half-plane above the graph y = arccose"^, x > and 
y = 0, X < 0. Notice that H = Hi, where Hi was defined in the proof of 
Theorem 1. 

Thus g{A) — h{A) -^ const where the constant is given by 

TT Jo t 7C Jo 2 t 

The integral is convergent because Vo{t) = 0{\/i), t — > and Vo(t) = 
0(e~*), t — >■ 00. It remains to find the asymptotic behavior of g{A) as 
A -^ 00. We have 

Combining these results we obtain 

B = h{A) = A + -log A+-J^ ^^H{t)--X[2,oo){t))j + o{l) 

= A + -logA + c + o(l), 
where c is the constant from (fT^ . This proves Theorem 2. 

Proof of Theorem 5. We choose 

B = log m — log log m + log 4, (16) 

so that L ~ (log?7i)/(2m) in view of (0)- Consider the conformal map of 
^B,m onto the half-strip U = {z : ^z & {0, {m + 1)11)} by a function /i such 
that /i(0) = 0, /i((m + l)7r) = (m + l)7r and /i(oo) = 00. Let B' = fi{B). 
The function fi{z/B) tends to the identity as m — ;> 00, so 

B' r^ B, B ^ 00. 

Let /2 be the (elementary) conformal map of n onto the first quadrant nor- 
malized by f2{B') = 0, /2(7r(m + !)) = ! and f2{oo) = 00. Then 

" ■= ^^^°) = ! + exp(fv(^t!)) = ^'/^^"^ + ^^ ^ ^/^^ ^ (logm)/(2m). 

The authors thank Doron Lubinsky, Misha Sodin and Andrei Gabrielov 
for their help and useful comments. 
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